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Abstract: We present an analysis of supersymmetric left-right symmetric models with 
Higgs fields lying in the adjoint representation of S'[/(2 )r. These models feature a doubly- 
charged Higgs boson which gets its mass only at the loop level. We present, for the 
first time, a complete one-loop calculation of this mass and show that contributions that 
have been neglected so far can shift it by a few hundreds of GeV. We then combine this 
observation with LHC bounds deduced from extra charged gauge boson and doubly-charged 
Higgs boson searches. In particular, we point out that existing limits get substantially 
modified by the presence of singly-charged Higgs bosons that are also predicted in these 
models. In addition, we address constraints stemming from vacuum stability considerations 
and show how the considered class of models could be ruled out at the next LHC run in 
the absence of any signal. 
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1 Introduction 

The Large Hadron Collider (LHC) at CERN has been built to explore the TeV energy scale 
and pin down the mechanism responsible for the breaking of the electroweak symmetry. 
With this respect, the discovery of a Higgs boson consistent with the Standard Model 
of particle physics [1, 2] can be seen as a first success of the LHC physics programme. 
This also consists of the first observation of a particle intrinsically unstable with respect 
to quantum corrections. One can therefore expect either an unnatural fine-tuning or a 
stabilization arising from physics beyond the Standard Model that would emerge at a scale 
reachable at present and future colliders. Although we do not know which theory could 
be the most suitable, a lot of efforts have been put, during the last decades, in the study 
of a plethora of new physics models. The Minimal Supersymmetric Standard Model [3, 4] 
(MSSM) is one of the most popular of these and it provides, additionally to curing the 
infamous hierarchy problem, a solution for other conceptual issues plaguing the Standard 
Model. 

There is currently no evidence for supersymmetry and limits on the superpartners of 
the Standard Model particles are hence pushed to higher and higher scales. Most results 
have however been derived either in the MSSM or within simplified models inspired by 
the latter. There are nevertheless large varieties of alternative supersymmetric realizations 
that could evade all current bounds and deserve to be investigated. For instance, the 
MSSM inherits some of the flaws of the Standard Model and extending the model gauge 
symmetry could provide a mechanism yielding the generation of the neutrino masses and 
an explanation for the strong CP problem. 

We consider in this work supersymmetric theories that exhibit a left-right symmetry [5- 
8] and that are based on an SU{3)c x SU{2)l x SU{2)r x U{1)b-l gauge group. They 


- 1 - 




provide a solution for the hierarchy problem, explain the smallness of the neutrino masses, 
can include a dark matter candidate, and they potentially solve the strong CP problem 
because of parity invariance at the fundamental level. Moreover, minimal supersymmetric 
setups like constrained versions of the MSSM cannot usually predict, at the tree level, a 
neutral Higgs boson mass compatible with the observations and require large higher-order 
contributions. This issue can be alleviated by extending the Higgs sector, like in left-right 
supersymmetric models where the spontaneous breaking of the left-right symmetry all the 
way down to the 17(1) electromagnetic group necessitates several Higgs fields [9-11]. 

Several left-right supersymmetric setups have been studied in the past, with different 
Higgs sectors and symmetry breaking details. We choose to focus on a class of mini¬ 
mal models featuring two SU{2)ji Higgs triplets that are responsible for the breaking of 
SU{2)ii X U{1)b-l to the electroweak group, two extra SU{2)i Higgs triplets that allow 
for parity invariance at higher scales and two SU{2)l x SU{2)ji bidoublets yielding both 
the breaking of the electroweak symmetry down to electromagnetism and the generation 
of acceptable fermion masses and mixings [12-14]. In addition, we include an extra singlet 
field to achieve a successful electroweak symmetry breaking in the supersymmetric limit [8] . 
With such a choice for the Higgs sector, the interactions of the lepton supermultiplets with 
the triplets are sufficient for generating neutrino masses through a see-saw mechanism [15], 
and the presence of the singlet solves the so-called //-problem of the MSSM as all bilinear 
terms of the superpotential can be generated dynamically. These models are currently 
probed at the LHC through traditional supersymmetry [16, 17] and extra gauge boson 
searches [18, 19], as well as analyses dedicated to the doubly-charged (Higgs and Higgsino) 
states induced by the SU{2)l and SU{2)b triplets [18, 19]. The observation of such doubly- 
charged particles would indeed provide valuable information about the symmetry breaking 
pattern of the underlying theory. 

The minimization of the scalar potential of the model is known to lead to a tree-level so¬ 
lution in which either i?-parity invariance [12, 13] or electric charge conservation [14] is lost. 
Including one-loop heavy Majorana neutrino contributions however allows for satisfactory 
solutions to the minimization equations [8, 20]. In this case, the right-handed sneutrino 
fields are prevented from getting non-vanishing vacuum expectation values so that i?-parity 
is conserved and no dangerous lepton number violating operators appear in the superpoten¬ 
tial. Additionally, these radiative corrections also protect the doubly-charged Higgs bosons 
from being tachyonic and breaking electromagnetism. This has been demonstrated both 
in pioneering calculations relying on the Coleman-Weinberg effective potential [8] and in 
a more recent calculation using a diagrammatic approach in the gaugeless limit where all 
contributions proportional to the gauge couplings are neglected [20]. 

In this work, we start with a brief description of the considered left-right supersym¬ 
metric models (section 2) and then show, in section 3, that the previously calculated 
corrections to the doubly-charged Higgs mass are not sufficient to render it compatible 
with the current experimental bounds [21]. We subsequently evaluate the complete one- 
loop result and demonstrate that the doubly-charged Higgs mass can be increased well 
above the present limits. Other recent LHC constraints are investigated in section 4 and 
we further study the minimum of the scalar potential at the one-loop level in section 5, 
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analyzing the requirements for ensuring this minimum to be global. Our conclusions about 
the phenomenological viability of the model are given in section 6 . 


2 The Model 


Many left-right supersymmetric setups have been developed in the past. We adopt the 
configuration introduced in ref. [ 8 ]. The Standard Model left- and right-handed quarks 
(leptons) are embedded, together with their scalar superpartners, in Ql and Qr {Ll and 
Lr) supermultiplets lying in the fundamental representation of the SU{2)r and SU{2)r 
gauge groups, respectively. The Higgs(ino) sector is constituted of two SU (2 )l triplets (Ai^ 
and A 2 L with a B—L charge Qb-l = T2), two SU{2)r triplets (Ai/^ and A 2 _r with Qb-l = 
=F2), two SU{2)l X SU{2)r bidoublets (d>i and ^2 not charged under the B—L symmetry) 
and one gauge singlet {S). As demonstrated in pioneering works, this configuration allows 
for a successful symmetry breaking of the SU(2)l x SU{2)r x U{1)b-l gauge symmetry 
down to the U{1) electromagnetic gauge group, it includes a see-saw mechanism yielding the 
generation of the neutrino masses and it provides acceptable fermion masses and mixings 
with respect to data [12-15]. 

Assuming a discrete Z 3 symmetry, the superpotential only contains terms that are 
trilinear in the fields, 

bb = QlUi^iQr + QLy2^2QR + LrU^^iLr + LLy2^2LR + A2lTl 
+ Lru^AirLr + XlSAil-A2l + ^rSAir-A2r + -|- A25'<h2-‘h2 ( 2 .I) 

-|- Ai25'<I>i •<I>2 -|- ) 


where we refer to ref. [22] for details on the underlying (understood) index structure. In 
those notations, the various coupling strengths have been embedded into 3x3 Yukawa 
matrices (t/q and with i = 1,2 and j = 1 , 2 , 3,4) and A parameters dictating the size 
of the Higgs self-interactions. The supersymmetry-breaking Lagrangian contains sfermion 
and gaugino mass terms as well as trilinear (T) interactions whose structure is given by 
their superpotential counterparts. 

The spontaneous symmetry-breaking mechanism proceeds in two steps. The SU{2)r x 
U{1)b-l group is first broken to the hypercharge group and the two SU{2)r Higgs triplets 
develop non-vanishing vacuum expectation values (vevs). In a second step, the electroweak 
symmetry is broken to electromagnetism and the bidoublet Higgs helds get non-zero vac¬ 
uum expectation values. The vacuum state is characterized by the vevs of all the neutral 
components of the Higgs fields. 


(Ail) = 




(Air) 


-^1 



(A2l) 

(A 2 r) 




A^r+ 

_^2R 

V2 
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( 5 ) = 


vs 
^/2 ’ 



where all vevs are taken real and positive (as allowed by suitable field redefinitions). Al¬ 
though the fields could also develop non-zero vacuum expectation values giving rise 
to W/Wn mixing, these are strongly constrained by kaon data so that they are ignored. 
Moreover, the phase of vs cannot in principle be rotated away but we omit it for sim¬ 
plicity. The remaining vevs can be further constrained. Both the ViR vevs are imposed 
to be large by the masses of the Wr and Zr gauge bosons. Furthermore, the small¬ 
ness of the Standard Model neutrino masses and electroweak precision data require that 
the vevs developed by the SU{2)r Higgs triplets vn are negligible^. Finally, the size 
of Vs is related to the effective ^-terms of the superpotential, and we consequently im¬ 
pose the hierarchy vs,vir,V 2 r 3> Vu,d S> uii,~U2l~0 . For further references, we dehne 
~ tanI3 r = V 2 r/vir, v‘^ + v'^ = and tan/3 = Vujvd and we additionally 
assume Ai = A 2 = 0, analogously to refs. [22, 23]. 

The vacuum configuration of eq. (2.2) does not necessarily yield a minimum of the 
scalar potential [12]. In fact, the true global minimum is in general charge-breaking and 
the and fields develop a vev 


/ 0 

(Aiij)cB=(^ Q j for z = 1,2, (2.3) 

so that the H-term contributions to the scalar potential are minimized. We will further 
address this issue in section 3. 

In order to design viable left-right supersymmetric benchmark scenarios, we have im¬ 
plemented the model described above in the SARAH package [24-28] to automatically gen¬ 
erate a numerical code based on the SPheno programme [29, 30]^. Physical spectra are 
then calculated in several steps. One begins with the derivation of the soft-supersymmetry 
breaking masses from the tadpole equations at the tree level, which allows for the calcu¬ 
lation of a tree-level mass spectrum. Next, one-loop corrections are evaluated in the DR 
scheme and a one-loop-accurate mass spectrum is computed. As explained in the next 
section, we have modified the automated SARAH-SPheno procedure due to special features 
associated with the mass of the doubly-charged Higgs boson. Several analytical cross checks 
with the model implementation in FeynRules [31, 32] have been performed. For the evalu¬ 
ation of the LHC bounds performed in section 4, the SSP package [33, 34] has been used to 
scan randomly the parameter space of the model. Moreover, we have estimated Wr signal 
cross sections by means of the MadGraph5_aMC@NL0 programme [35] after having linked the 
UFO library [36] obtained from the SARAH implementation. 

^Another acceptable choice implies that the j/3 couplings are tiny and the j/f ones are large. This is 
however very unnatural with the left-right symmetry. 

^The code can be obtained from the authors upon request. 
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Figure 1. The squared mass eigenvalue of the lightest doubly-charged Higgs state for 

the desired vacuum configuration of eq. (2.2) (left panel, blue) and for the charge-breaking one of 
eq. (2.3) (left panel, purple) as a function of tan/?i^. We also zoom into the region with tan/3i^ « 1 
(right panel) and present results for the tan /? = 1 (blue) and tan /3 = 50 (red) cases. Moreover, the 
condition > 2|u|^ — v'lp\ that is discussed in the text is illustrated by vertical lines. The 

relevant model parameters have been set to vr = 5.5 TeV, vs = 10 TeV, Xr = 0.5, Xr = —0.5, 
Ai 2 = —0.02 and T\^ = 0. 


3 The doubly-charged Higgs mass 

In the class of supersymmetric left-right symmetric models under consideration, we have 
four doubly-charged Higgs bosons, two of them being related to the SU{2)r sector and two 
of them being related to the SU{2)r one. While at the tree level, these fields do not mix, 
a small mixing is induced at the loop level. It however does not lead to any observable 
effect so that the following discussion is solely focusing on the SU{2)r sector. The mixing 
effects are nevertheless included in our numerical results. 


3.1 Tree-level results 

The squared mass matrix of the doubly-charged SU{2)r Higgs bosons reads, at the tree 
level and in the (A)”^, A^^*) basis, 

~ tan/3ij F++ F++ 

^ y F++ -D++ - cot fin F++ 



In our conventions, the H-term and T-term contributions are respectively given by 


D++ = 


4 ^ 

2 L 


vl-vl + - vIr) 


JP 4 , 2 

F++ = -y'^^R'^'^R - 


ArAi2 


(3.2) 


where qr stands for the SU{2)r gauge coupling. First, the form of this squared mass 
matrix implies that one of the physical states is massless in the gaugeless limit, i.e., when 
all gauge couplings are neglected. Next, once the contributions proportional to the gauge 
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Figure 2. Generic one-loop diagrams contributing to the doubly-charged Higgs self-energy. 

couplings are taken into account, the condition 

bu “(3-3) 

must be satisfied in order to render both doubly-charged Higgs squared mass eigenvalues 
positive [14] . Equivalently, this condition can be casted as tan /3ij — 1 <C 1 since we have 
imposed, for phenomenological reasons, that vr ^ Vu^Vd (see section 2)^. 

The condition of eq. (3.3) turns out to be extremely fine-tuned, as illustrated in figure 1 
where we depict the dependence of the smallest mass eigenvalue on tan /3r. In the left 

panel of the figure, we present numerical predictions for the desired vacuum configuration 
of eq. (2.2) (blue) as well as for the charge-breaking one of eq. (2.3) (purple). In contrast 
to the latter case where the squared mass is generally positive, there exists only a small 
region of the parameter space, featuring tan/?/j « 1, where the vacuum structure does not 
break electric charge invariance. In the right panel of hgure 1, we zoom into this region and 
study scenarios for which tan /Sr ~ 1. In the case where tan /? = 1 (blue), or in other words 
when — v‘^\ = 0, the condition of eq. (3.3) imposes that the lightest doubly-charged 
Higgs boson is either tachyonic or massless. In contrast, when tan (3 increases, a small 
window where the squared mass of the lightest doubly-charged Higgs boson is positive 
opens (red, for the example of tan/3 = 50). However, this conhguration implies that one 
of the pseudoscalar Higgs bosons gets tachyonic [37]. 

As a consequence, it is impossible to construct, at the tree level, left-right supersym¬ 
metric setups that are phenomenologically viable. The only alternative option would be 
to consider solutions of the potential minimization equations featuring i?-parity breaking 
vevs for the right-handed sneutrinos [12, 37]. 

3.2 Results at the one-loop level 

The problems mentioned above are a tree-level artefact and can be solved by considering 
one-loop contributions. Attempts including quantum corrections derived either from the 
Coleman-Weinberg potential [8, 38] or within a Feynman diagrammatic approach [20] have 
shown that a charge-conserving minimum of the scalar potential could be obtained. A 
complete one-loop calculation however exhibits logarithmic terms with a negative argument 
resulting from contributions depending on the self-interactions of the doubly-charged Higgs 
bosons, and the subsequent imaginary part of the one-loop corrected potential is a sign 

®The special case where tan Pr — 1 corresponds to a saddle point and not to a minimum of the potential. 
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Figure 3. Comparison of our calculation of the one-loop corrections from j/f, evaluated in the 
gaugeless limit and for one generation of right-handed neutrinos (red dashed), to the results of 
ref. [20] (grey). The contour lines represent isomass lines in GeV^. We present results 

in the (wij, logiQ(m|_^/GeV^)) plane with yf = 0.4 (left panel) and in the (uij,yf) plane with 
= 2 • 10® GeV^ (right panel). In both cases, the remaining relevant parameters are fixed to 
Xr = 0.4, tan/3fl = 1.02 and vs = 10 TeV and the white areas are regions where one of the right 
sneutrinos is tachyonic. 

of a dangerous unstable vacuum (see section 5). More precisely, these earlier works have 
argued that contributions stemming from the Majorana neutrino Yukawa coupling were 
sufficient to render the doubly-charged Higgs squared mass positive [20] and to make the 
vacuum configuration of eq. (2.2) a deeper minimum than the one of eq. (2.3) [8]. However, 
in this last case, the right slepton soft squared masses must be negative, which further 
implies large one-loop gaugino contributions to the slepton mass eigenvalues to avoid charge 
and i?-parity breaking. 

In this work, we have extended the Feynman-diagram-based calculations of ref. [20], 
the latter having been achieved in the gaugeless limit, for one generation of right-handed 
neutrinos and considered only diagrams depending on y^. We have instead computed all 
one-loop contributions to the doubly-charged Higgs boson mass, the associated Feynman 
diagrams being presented in figure 2. In figure 3, we start by numerically comparing our 
results, evaluated in the appropriate limit, to those of ref. [20] and demonstrate that a good 
agreement has been found^. In the parameter space regions probed on the figure, the soft 
slepton mass and the effective supersymmetric Higgs mass parameter = Xrvs / 

are not too large, which leads to a doubly-charged Higgs boson that is still tachyonic. The 
loop corrections are even in this case negative and thus counterproductive to restore a 
vacuum state that conserves the electric charge. Only for large values of Xrvs (not shown 

^The authors of ref. [20] have confirmed a typo in the first line of their eq. (138) in which the should 
be read outside of the logarithm. 
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on the figure), a positive contribution could emerge. This requirement however also yields 
a significant CP-splitting of the right sneutrinos, so that large values for are as well 
required to prevent them from being tachyonic and breaking P-parity. Assuming that all 
sfermions have similar masses, realistic setups of the class of models under study are thus 
unlikely to be observable from standard sfermion searches at the LHC. 

As shown in the rest of this section, the situation improves when the complete one-loop 
result is considered. Technically, our predictions rely on the SARAH-SPheno procedure for 
generating the particle spectrum. In the usual running mode of both programmes, the 
tree-level DR-masses are used for the particles running into the loop diagrams. In the left- 
right supersymmetric models under consideration, this leads to problematic log(m^±±/Q^) 
terms as the lightest doubly-charged Higgs boson is tachyonic at the tree level. To avoid 
this issue, we have modified the SARAH-SPheno procedure in a way that is inspired by the 
on-shell scheme, and inserted one-loop corrected masses into the loops. Iteratively, this 
consists of: 

• hrst iteration: using = |[m|^±±]tree-ievei|; 

• ^\k] = until < e < 1 . 

We have fixed e = 10“^ so that four iterations are generally necessary for any benchmark 
scenario. Moreover, the initial value is in principle arbitrary so that we could have 
chosen = 0. We have verified the invariance of the predictions with respect to this 
choice. Although our prescription breaks gauge invariance at the two-loop level, the associ¬ 
ated effects are expected to be significantly smaller than the genuine two-loop contributions 
and thus under good control. Here and in the following sections mjj±± denotes the mass 
of the lightest doubly-charged Higgs boson. 

In figure 4, we show predictions for the mass obtained when using our complete 
one-loop calculation. The results are presented in logj^olui'i^/GeV^)) planes, first (left 
panel of the figure) when all the other model parameters are fixed as in the left panel 
of figure 3. Although both the tree-level and the approximate one-loop predictions yield 
tachyonic doubly-charged Higgs bosons, the fully-one-loop-corrected mass is positive and of 
the order of a few hundreds of GeV. Moreover, the dependence on is explained by the 
2 / 4 -dependent contributions. In the right panel of figure 4, we increase the value of Xji and 
decrease the one of so that the dependence on almost vanishes and the one-loop- 
corrected mass reaches values beyond 500 GeV. The most important contributions 

to the full one-loop result consist of positive corrections arising from Wr/ singly-charged 
Higgs boson and doubly-chargino/neutralino loop-diagrams, as well as from negatively 
contributing yf-induced loop-diagrams, singly-charginos and doubly-charged Higgs bosons. 

The parametric dependence of the loop contributions is highly non-trivial because of 
the large number of different contributing sectors. However, some generic features can be 
extracted using the information given in figure 5, where we add successively different loop 
contributions for an exemplary benchmark scenario. We start by adding to the tree-level 
results the yf-dependent contributions which have already been discussed: they become 
negative for large values of yf and 
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Figure 4. Contours of mH±± [in GeV] including the complete one-loop calculation in the 
(u_R, logig(TO|^/GeV^)) plane for a scenario featuring the same setup as the left panel of figure 3 
(left), as well as for a scenario featuring a smaller yf = 0.1 and a larger Xu = 0.9 value (right). 


In order to describe the additional contributions that will be further added, we recall 
that the structure of the relevant chargino and neutralino couplings reads 


-pL 



1 , 

■pi? 

xt+PgH-- 



> (3.4) 

pL 

xtx+H^- 




pi? 

xtXjH-- 





where we have parameterised the vertices as F = In our notations, U^, 

and Z are the mixing matrices of the neutralinos, charginos and doubly-charged Higgs 
bosons respectively, whose elements are written as Ui^x and Zi^x- These matrix elements 
hence denote the X field component of the mass eigenstate. We also introduce a set of 
vector and Higgs boson interactions whose strengths are given by 


^H--W+W+ 

^H--H+W+ 


iV2gl{v,nZ-^ +v,nZ-^ ) 






(3.5) 


-(z 


— ~^9r (2 ( V2rZ . i Z 


. . + —, . + Z, .++,* -|- virZ, „Z „Z. — 


^PlR jPlR ^p2R 


V2rZ. .— -h VirZ 

R 


2, A, 


+ + )* 


the momentum and metric dependence being omitted from the Feynman rules. 
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Figure 5. Dependence of on the S'C/(2)fl-breaking scale vr. The lines correspond to to 

the tree level prediction (grey solid) and add then successively different contributions: tree-level 
-I- yf-dependent (s)lepton and (s)neutrino contributions (blue dotted), + neutral/doubly charged 
Higgses and neutral gauge bosons (blue dashed), -I- chargino contributions (blue solid), -|- neutralinos 
and doubly-charged Higgsinos (black dashed). The full result (black solid) finally contains also 
Wr/H^ contributions. The employed benchmark scenario is defined by Xr = 0.4, tan/J^j = 1.02, 
"Ilh = 2 • 10® GeV, j /4 = 0.25 and features one generation of right-handed neutrinos. 


Beside the contributions, a set of self-energy diagrams also contribute negatively. 
These involve either the doubly-charged Higgs boson and a neutral vector or Higgs boson 
or two singly-charged charginos. The sign of the latter, where one would have naively 
expected a positive contribution, can be understood from the large mixing between the 
gaugino and the singly-charged Higgs bosons of the SU{2)r sector. The contributions 
to the diagonal entries of the doubly-charged Higgs mass matrix are proportional to 
(|r{'I .1 _P-|-|r{^i .1 „_whereas the contributions to its off-diagonal entries are 

' ' xJXjHf^ ' ^ 

proportional to (T'^', „_)(r(^, „_)*7Ti-+m-+. The masses of the relevant charginos 

are however proportional to vr and the wino soft supersymmetry-breaking mass so that 
they are in general much larger than mj{±±. The contributions to the off-diagonal entries 
of the doubly-charged Higgs mass matrix turn therefore out to be much larger than those 
to the diagonal ones. After diagonalisation, this yields a negative contribution to mR±±. 

The main positive contribution to the doubly-charged Higgs mass are loop diagrams 
containing a neutralino and a doubly-charged higgsino. The difference with the singly- 
chargino case depicted above stems from the singlino component of the neutralinos. In 
the respective entries of the mass matrix, the product (r^{+ _)(r'^{+ _)* involves 

xVx'jH^^ XVX)H^ 

Ai^-dependent terms whose sign is different from the one of the gauge contributions. These 
terms will hence dominate for large values of Xr. Additionally, the Wr/H~^ loop diagrams 
also contribute positively to the diagonal entries of the mass matrix. 

Eventually, the negative growth of the tree-level contribution to the mass (grey 
line in figure 5) with increasing vr cannot be compensated anymore so that the gets 
again tachyonic for large enough vr. This scale turns to be vr ~ 14 TeV for the example 
of figure 5. 

Both LHC collaborations have set bounds on the mass of the doubly-charged Higgs 
boson. These limits however strongly depend on the final state in which the doubly-charged 
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Higgs boson decays into [21, 39]. In this way, mfj±± is constrained to be larger than 
204 GeV, 459 GeV, 396 GeV and 444 GeV for rr, ///r, e/r or ee final states, respectively. 
As shown in figure 4, there is usually a large range of vr values, while keeping all other 
parameters fixed, for which these LHG constraints are satisfied. However, vr also sets the 
mass scale of the Wr boson so that one can combine these constraints with LHG limits on 
the Wr mass (see section 4). 

4 LHC bounds 

The model parameters entering in the above computations can also be experimentally 
constrained by searching for the rest of the particle spectrum. We focus here on the 
electroweak part, i.e., extra gauge bosons and right-handed neutrinos, updating the results 
first given in the report of the 2013 Les Houches workshop [23]. 

In the left-right supersymmetric setup under study in this work, the extra Zr boson 
is always heavier than its Wr charged counterpart. We hence focus on the latter since it is 
more constraining. In general, a Wr boson can decay into pairs of fermions as the Standard 
Model W boson, into pairs of Standard Model bosons, into pairs of sfermions, into a 
chargino and a neutralino final state, and into a charged lepton and a right-handed neutrino. 
Therefore, one can experimentally look for it in several search channels, each setting its own 
bounds. The seemingly most stringent one is typically set by investigating the signature 
of a charged lepton and missing energy, assuming that the new Wr boson decays into a 
lepton and a low-mass neutrino which escapes detection. However, in a left-right scenario 
this bound does not apply. When neutrino data is explained by a seesaw mechanism, the 
decay of the Wr into a lepton and a low-mass neutrino is generally suppressed by the small 
neutrino mixing, while the Wr —)• Iur mode will typically give rise to more complicated 
final states. Before focusing on the latter channel, more involved due to the unknown 
right-handed neutrino mass, we review the simpler searches with hadronic two-body final 
states, i.e., Wr —>• jj and Wr —)■ tb. In the first case, we can reinterpret the CMS inclusive 
dijet search of ref. [40]. In the second case, a direct comparison to the CMS search for 
Wr bosons in the tb channel of ref. [41] is in order. We do not consider the ATLAS 
counterparts because they are less sensitive. The above searches do apply to left-right 
models since the coupling strength of the Wr boson to a pair of quarks is equal to that of 
the Standard Model W boson, given that qr = gi holds. Their reinterpretation is done by 
comparing the pp —?■ Wr —)• jj{tb) cross sections evaluated with the MadGraph5_aMC@NLD 
programme [35] to the excluded cross sections as a function of the Wr mass. More precisely, 
our predictions have been evaluated by convoluting leading-order squared matrix elements 
with the CTEQ6L1 set of parton densities [42] and include a AT-factor of 1.23 and 1.2 for 
PP jj and pp —)• tb, respectively. 

In figure 6 the cross sections excluded by CMS are compared to those computed in 
our model. In the right panel of the figure, we observe that considering the full model 
spectrum decreases the constraint by typically 50 — 100 GeV, due to the appearance of 
extra decay modes for the Wr that are not considered in the simplified models adopted 
in the experimental analyses. These extra channels get the Wr decay width proportion- 
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Mw^ [GeV] Mw^ [GeV] 

Figure 6. Bounds on the boson obtained from the CMS dijet search (left) and analysis of 
tb leptonic events {£ = e, /i, r) (right). The CMS exclusion curves are taken from ref. [40] and 
ref. [41], respectively. In the latter search, CMS puts a direct bound of Mwr > 2.05 TeV, while no 
such bound exists in the former CMS analysis. 

ally larger, slightly suppressing the branching ratios into the hnal states that are searched 
for. The spread in the model prediction reflects the scan over the parameter space. The 
dijet hnal state, shown in the left panel of the hgure, is the one setting the tightest 
bound, Mp/jj > 2.3 TeV. The corresponding bound on the right-handed vev is given by 
vr > 4.9 TeV. A value of the Wr boson mass between 1.8 and 2.0 TeV, which could be 
allowed by the dijet search depending on the specihc choice of the input parameters, is 
fully excluded when the tb search is also considered. 

We move now to the constraints in the Wr —)• Ivr process. It is usually assumed in 
left-right symmetric models that the right-handed neutrinos can only decay via the Wr 
boson, which can be either on- or off-shell, with the subsequent decay of the latter into two 
jets in two thirds of the cases. Therefore, footprints of the left-right models are searched 
for in the pp — )• Ujj {I = e, p) hnal state. This is used by experimental collaborations 
to set strict bounds on a combination of masses of and Wr, as, e.g., in ref. [43]. As 
shown in ref. [23, 44], this simplihed assumption can sometimes be too restrictive. This 
is particularly true in the left-right supersymmetric models under consideration, where 
new two-body decay channels for the ur can be present, such as ur —)• and decays 

into a neutral (charged) slepton and a chargino (neutralino). Furthermore, several other 
decay modes for the Wr boson exist, thereby providing new three-body decays for the 
right-handed neutrinos. This is pictorially shown in hgure 7. 

All these new decay modes for the right-handed neutrinos have the net effect of reducing 
its branching ratio into the search channel of interest, i.e., vr — )• £jj. Two main effects are 
visible. The right-handed neutrino prefers to decay into the lightest charged Higgs, if this 
channel is kinematically viable, with the subsequent H'^ —)• tb decay. However, this decay 
is possible only due to the (small) right-handed Higgs triplet component of the lightest 
charged Higgs. Therefore it is dominant only when the Wr needs to be far off-shell for the 
three-body decays to occur, otherwise the latter will dominate, as happens for large masses 
of the right-handed neutrino. We observe the second effect on the hgure. The presence 


- 12 - 














400 < MH±/GeV < 500 


rrr^^ other I 


500 1000 1500 

[GeV] 


2000 


Figure 7. Right-handed neutrino branching ratios for a charged-Higgs mass of 400 GeV < M^± < 
500 GeV. ‘Other’ refers to all other three-body decays via an off-shell Wr boson. 
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Figure 8. Allowed (green) and excluded (red) points from the CMS search for Wr and i^r in the 
eejj channel for a lightest right-handed neutrino of electron flavour (of mass for 400 < 

MR±/GeV < 500, in the (m^^^,Mn/^) plane (left) and (j/|^[l, 1],u/j) plane (right). The CMS 
exclusion curve (in black in the left panel of the figure) is taken from ref. [43]. The red shaded area 
represents the dijet bound. 


of extra decay modes for the Wr (mainly charginos and neutralinos) further suppresses 
the searched channel, which in the experimental approximation has a constant branching 
ratio of 67%. However, the branching ratio —)• £^H^) depends on the right-handed 

neutrino mass as well as on the singly-charged Higgs mass. It can get up to 0(80%) for 
very light singly-charged Higgs bosons and right-handed neutrinos and conversely heavy 
Wr bosons. 

As done previously for the hadronic two-body decays of the Wr boson, we now reinter¬ 
pret the CMS bounds on the Wr —)■ Irr —>• iijj channel of ref. [43]. Figures 8-9 display the 
allowed points (green dots) as well as the excluded points (red dots) when the considered 
right-handed neutrino is of electron and of muon flavour, respectively, in the (mi/^,Mu/^) 
plane on the left panel of the figure, and in the {y\,VR) plane on its right panel. 

We have found that the bounds set by CMS are generally too strong to be naively 
imposed in our left-right supersymmetric setup, and quite large areas within the excluded 
regions survive, especially for an electron right-handed neutrino. Investigating also figure 7, 
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Figure 9. Same as in figure 8, for a lightest right-handed neutrino of a muon flavour (of mass 

we identify three regions where exist parameter configurations which are allowed. First, 
for very light right-handed neutrinos, its decay products are too few energetic and appear 
below threshold. This consists of the lower region, below the black curve, in figures 8-9. 
Next, the region at high masses for both the Wr and the vr, that contains allowed points 
because of extra two-body decays of the Wr that suppress the search channels. Last, the 
bulk area in figure 8 is now allowed, contrary to the CMS results, because of the suppression 
of the associated cross section due to the vr —?■ decay channel which is dominant 

for moderate ur masses when such channel is kinematically open. Such an area does not 
appear for a muonic right-handed neutrino because of the much stronger CMS bounds. 
In fact, a fluctuation of about 3cr has been found in the eejj CMS search channel, for 
~ 2.0 TeV, which our model can fit, as noted in refs. [45, 46] for similar cases. 

5 Analysis of the vacuum stability 

In section 3 we have shown that loop corrections have the power to switch the SU{2)r x 
U{1)y breaking saddle point into a minimum and to give a positive squared mass for 
. Having now a local minimum with a proper breaking of SU{2)l x U{1)y does not 
necessarily imply that this is also the global minimum of the theory. In particular, it is not 
obvious that it is deeper lying than the charge-breaking minimum of eq. (2.3) after taking 
into account consistently all loop effects. For further references, we denote in the following 
the vacuum configuration with the desired symmetry breaking pattern as ‘DSB’ and the 
charge-breaking vacuum as ‘CB’. 

The tadpole conditions cannot be solved analytically at the one-loop level so that one 
has to rely on numerical methods. However, the main features of the vacuum configura¬ 
tion can already be understood by considering a simplified model where we focus on the 
SU{2)r X U{1)b-l gauge sector, including only the two SU{2)r Higgs triplets Air and 
A 2 R as well as one generation of right-handed (s)leptons. The corresponding superpotential 
is given by 

^ = yA Lr Air eLR + XrS Tr(Ai/j A 2 _r) + - A 5 5^ , (5.1) 


400 ^MH±/GeV < 500 
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where we have written down the full SU{2) matrix structure for the first term using 
Lr = {^rA^r}^ as well as the e tensor defined by ei 2 = —621 = 1. The tree-level scalar 
potential Vq reads 

Vq = Vf + Vd + 14 oft ( 5 . 2 ) 

with 


Vf 


Vd 


Koft 


\l |5|2 (Tr(Al^ Ai^) + Tr{Al^ A2 r)) + | Ar TV(Air A 2 r) + A 5 5' 


+ 2/4 Ar ( S* Lr A 2 ^ e Lr + h.c. 


+ ( 2/4 


Y ^ 1 ^]) + Tr(A^^[Ti, A2 r]) - L^jTi Lr) 

i=l 

+ ^ ('iv(At^ A2 r) - TV(a1^ Air) + ^lI^Lr) 

”T'AiflTr(AiR A\j^) + mi^^Tr(A2R A|^) -h m|^L^LR 

+ (Ta^ 5 TV(Air A 2 r) + + Ti Air e Lr + h.c.) . 


( 5 . 3 ) 


In these expressions, r* = ^cjj, where dj are the Pauli matrices. For simplicity we set 
the trilinear soft supersymmetry-breaking couplings to zero as they do not change the 
qualitative features that we are aiming to describe. 

We start by removing the soft supersymmetry-breaking masses and 

from the equations by requiring that the tadpole equations for the desired charge-conserving 
case of eq. ( 2 . 2 ) are solved, 


m 

dX 


DSB:(A0^> = ^,(S> = );| 


= 0 


‘x /2 


for X = {A?^,A0^,5}, 


( 5 . 4 ) 


so that the expression of Psoft is optimized for the DSB case. However, as argued above, 
the global minimum of Vq consists of a configuration where the vevs of the triplet fields 
are aligned along the ri direction. Taking the effect of the soft supersymmetry-breaking 
masses into account, the magnitudes of vir (with i = 1 , 2 ) get slightly modified by factors 
ai which are close to one, so that we can relate the vevs derived in the DSB case to those 
derived in the CB case. 


(A?r)|cb = )|cB = “ 1 “^ and (A^^)|cb = (A+^)|cb = “2^ ■ ( 5 . 5 ) 


Consequently, we trade the tadpole equations for the doubly-charged fields with 

^kblcB _ Q 
dui 


( 5 . 6 ) 


For the evaluation of the mass spectrum at the DSB and the CB minima, we split all 
complex scalar fields into their scalar and pseudoscalar components as for the CB minimum, 
as electromagnetism is eventually broken. Taking into account all possible vevs, we rewrite 
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every scalar field X as X = -^{vx + (p^+i4'x)- After rotating out the unphysical would-be 
Goldstone bosons in each configuration, we compute a 15 x 15 scalar mass matrix 


(^2)DSB/CB 


9 Vo 

d(l)id(j)j 


DSB/CB 


(5.7) 


The fermionic part of the spectrum can be evaluated for the SU{2)r {U{1)b-l) gauginos 
{B), Higgsino triplets as well as the the singlet fermion S and the lepton doublet 
Lb from the Lagrangian terms 

= - As 5 5 5 - (5 Tr(AiH X2r) + S Tr(Ai^ As/j) + 5 Tr(Ai^ As/j)) 

—^ Ai/j e Lr + Lr Air e Lr + Aijj e Lr^ 

3 

- V29R'^WR^i(Tr{A\j^[Ti, Air]) + TT{Alj^[Ti, A2 r]) - V^nlR^ 

2 = 1 

- V2gBLB (Tr(A^j^ A2r) - Tr(A|^ Ai/j) Lr^ 

1 ^ ^ _ _ X 

“ 2 ^Wr Y1 - 2 ^b BB + h.c.. 

2=1 

As an example for the differences between both vacuum configurations, the lepton masses 
are given by 

w-i/fllosB = V^viRy^, 
iRurIcb = «! viRui, 


W-eflloSB = 0 , 
IRejilcB = «! ViRH^ 


(5.9) 


Finally, the masses of the vector bosons are derived from the non-derivative part of the 
gauge-invariant kinetic terms of the Higgs bosons, 


=Tr ((5K E - 9bl B^ A\^){gR ^ [r^ Ai^] - gBL B, A,r) 

a=l 6=1 

3 3 

+ {qr E ^r" + 9bl B^ Alj^){gR E Wr^^ [r^ A2ij] -f gBL B^ A2 r)^ . 

a=l 6=1 

(5.10) 


This gives three heavy states of masses of 0{vr) and one massless state for each vacuum 
structure. In the DSB case, the hypercharge symmetry group remains unbroken by the 
triplet vevs and the associated boson is thus massless. In the CB case, the vevs are aligned 
along the ti direction so that this generator remains unbroken ([n, Aj^] = 0) and the Wr^i 
boson turns out to be massless. This is a consequence of the fact that a Higgs field in the 
adjoint representation of an SU{N) group cannot break the rank of this group. 

We now move on with the analysis of the one-loop effective potential 


blL,eff. = Ho -I- Vcw , 


(5.11) 
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tan i3r vr [GeV] 

Figure 10. Differences of the scalar potential depths between the DSB and the CB cases at the 
one-loop level. For negative (positive) values of this difference, the DSB (CB) minimum is the 
global one. In the left panel of the figure, we represent this difference at the tree (solid grey line) 
and at the one-loop (solid black line) level as a function of tan/3i^ for vr = 5.5 TeV, Xr = 0.4, 
2/4 = 0.25, = 2 • 10® GeV^ and vs = 10 TeV. In addition we split the results in terms of the 

various contributions: slepton and lepton fields (blue dotted line), Higgs fields (blue dashed line), 
gauginos/Higgsinos (red dotted line) and vector bosons (red dashed line). In the right panel of the 
figure, we depict the total difference for tan/^Tj = 1.02, i/f = 0 and Xr = 0.3 (solid black line), 0.4 
(green dotted line), 0.5 (blue dashed line) and 0.6 (red dot-dashed line). 


where Vcw is the Coleman-Weinberg potential. It reads, in the DR scheme and using 
Landau gauge, 


Vcw = ^ 

n 


647r2 




(5.12) 


where n runs over all real scalar fields, Weyl fermions and gauge bosons. We denote by Sn 
and rUn the respective spin and mass of the field and have also set the renormalization 
scale Q to 1 TeV. We show on the left panel of figure 10 the difference 


AV = V 


DSB 

lL,efr. 


- V, 


CB 


lL,efr. 


(5.13) 


at the tree-level, at the one-loop level as well as after breaking down the result for the dif¬ 
ferent contributions to the potential. For the chosen set of parameters, the global minimum 
is of the DSB kind up for tan/3i/ values in the [0.97, 1.03] range. As illustrated on this 
figure, the observed behaviour is a consequence of the interplay between the fermionic and 
the bosonic contributions in the Higgs and gauge sectors. In particular, the DSB vacuum 
can be the global minimum only due to the fermionic contributions. This may seem to 
contradict the results of the previous section where the charged Higgs and Wr diagrams 
are very important for getting a non-tachyonic doubly-charged Higgs. We however recall 
that we are focusing here on the differences between the minima and not on the absolute 
contributions. 

We have shown in section 3 that a sizeable Xr is needed to get positive. On the 

right panel of figure 10, we study the dependence of AV on vr for different values of Xr. 
We observe that increasing Xr also increases the relative depth of the DSB vacuum with 
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respect to the charge-breaking one. We have finally checked that the features discussed 
so far do not depend on the choice of Q by varying it up to 2vii and evolving the model 
parameters according to the renormalization group equations. 

We now turn to the study of the full model and allow for additional vevs. We have used 
SARAH to produce a Vevacious code [48] for the left-right supersymmetric models under 
study. The Vevacious programme starts by evaluating the tree-level scalar potential for 
a given spectrum and finds all extrema in terms of all scalar particle vevs, given together 
with the corresponding value of the potential. In a second step, the nearby extrema are 
found and evaluated using a full one-loop potential. In principle one should allow for all 
scalars to obtain a vev which however in models like the one considered in this work would 
take an enormous amount of CPU time well beyond a year for a single point of parameter 
space. Therefore, we restrict ourselves to vacuum configurations in which only the 4 * 12 ; 
‘^iR 2R ^ neutral fields and the A]”^ and doubly-charged fields could develop 

non-vanishing vacuum expectation values. We further allow for one generation of right 
sneutrinos to receive a vev in order to cover the possibility of a charge-conserving but 
i?-parity violating vacuum. 

Under these assumptions, in the parameter regions where an i2-parity conserving spec¬ 
trum can be found at the one-loop level, the global minimum is found to be always either 
of the desired kind of eq. ( 2 . 2 ), such that 


(ck?)= Vd/V2, 

{^ir) = vir/V2, 


(4>°)= Vu/V2, {S)= vs/yj2, 


(A 


\R 


= 0 . 


{v<^)= 0 , 

(5.14) 

(A^k)= 

(5.15) 

= 0 , 

(5.16) 


(5.17) 

parameter space 

where 


or the charge-breaking nature of eq. (2.3) with 

(<!>?) = (<kO)=0, (5)=us/V2, 

(A?r) = (A") ~ uik/2 , (A^^) = (A++) ~ V2rI2 . 


the desired vacuum conhguration corresponds to the global minimum of the scalar potential 
(green) and the ones where the charge-breaking configuration is preferred (red). We have 
found that regions relevant for LHC physics (he., regions where the mass of the doubly- 
charged Higgs is above the current limits) are not only exhibiting a viable local minimum 
but also a global viable one. Moreover, as for the 2 / 4 -dependent loop corrections to the 
mass, larger the and vr values are, the less favoured is the benchmark scenario under 
consideration. Furthermore, as illustrated in the right panel of the figure, a larger value of 
XrVs has the virtue of pulling the desired vacuum conhguration into a deeper minimum. 
However, as already above-mentioned, this increases the possibility of a tachyonic right 
sneutrino due to large negative T-term contributions proportional to the Xrvs product. 
This happens at the tree level in the parameter space above the black line. In the region 
indicated by light grey shading we consequently hnd i?-parity violating global minima, 
i.e. {D^) 7 ^ 0 , which still conserve electric charge but feature too large electroweak vevs of 
the order of one TeV. In the dark grey area we hnd minima which break both, i?-parity 
and electric charge. We stress that this i?-parity violation is solely due to the negative 
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Figure 11. Analysis of the vacuum stability. The results are presented in the plane with 

Xr = 0.3 and for one generation of right-handed neutrinos (left) and in the (j/f, Xr) plane with vr = 
5.5 TeV (right). The other model parameters have been fixed to tan Pr = 1.02, = 2 • 10® GeV^ 

and Vs = 10 TeV. The white contours depicted on the figures indicate isomass lines for the doubly- 
charged Higgs boson in GeV as obtained with our full one-loop-corrected calculation. The green 
regions correspond to setups where the vacuum configuration of eq. (2.2) is the global minimum of 
the scalar potential whereas in the red regions, the charge-breaking vacuum configuration consists 
of a deeper minimum. The blue lines separate the cases of a deeper DSB minimum and a deeper 
CB minimum in the context of the simplihed model discussed above. The parameter space above 
the black line features a tachyonic sneutrino at the tree level and scenarios in the grey regions 
consequently violate i?-parity while conserving (lighter grey) or breaking (darker grey) electric 
charge. 


F-term contribution which drives a sneutrino tachyonic and could be evaded by a larger 
soft slepton mass. It is therefore of completely different origin than the i?-parity violation 
suggested in early studies of this kind of models [12, 37]. Finally, considering more than one 
right-handed neutrino leads to a reduction of both the favoured region and the mass of the 
doubly-charged Higgs boson. The blue lines in figure 11 show the separation between the 
red and the green area as evaluated from the simplified model discussed in the beginning of 
this section. This simplified description can thus safely be used to get rough information on 
which of the two minima is the deeper one. It does however not cover the case of i?-parity 
violating global minima as one would need to extend the simplified setup with vevs for the 
sneutrino fields. 

Scenarios for which the global minimum is charge-breaking could nevertheless be viable 
in cases where the tunneling time from the local SU{2)l x U{1)y minimum to the global 
one is sufficiently large. The decay of a vacuum state into a deeper minimum can be 
described by the phase transition in which bubbles of true vacuum nucleate out of a false 
vacuum state. For a successful phase transition, the bubble has to be of a critical size and 
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(5.18) 


can, at zero temperature, be found by minimizing the Euclidean action [49] 

Se = J + V{4>)'^ , 

where cj) is a vector of all scalar fields and the coordinate corresponds to the imaginary 
time coordinate r = it. Equivalently, the equations of motion have to be solved so that 
6Se = 0, which eventually determines the optimal tunneling path in held space. For a 
hrst estimation of Se, the straight tunneling path between two vacua can be used, which 
corresponds to a reduction of the problem to a one-dimensional problem. The latter can 
be solved numerically to arbitrary precision by the so-called overshoot/undershoot method 
(see, e.g., ref. [50]). The full path deformation in all held dimensions is computationally 
extremely expensive and is hardly applicable to the dimensionality of the model studied in 
this work. The decay rate of the false vacuum per unit volume V is then given by [49] 

T/V = Ae-^^ . (5.19) 

In this expression, A is a quantity whose dimension is in the fourth power of energy 
and is related to the eigenvalues of a functional determinant [51]. This determinant is in 
practice usually estimated on dimensional grounds [50] , A ~ with M being the typical 
mass scale of the model. 

A hrst estimate using the direct tunneling path at zero temperature shows that all 
scenarios included in the excluded region of the hgure (red) are metastable which is due to 
the large separation in held space of the different vacua. 

It is however well known from the MSSM [52, 53] that allowing for additional vevs 
and/or including thermal effects can imply that metastable vacua are in fact unstable. A 
thorough investigation of this effect lies beyond the scope of this paper. 

Combining the results of the right panel of hgure 11 with the LHC constraint of 
vr ^ 4.9 TeV (see section 4) implies that all but a small strip of the parameter space turns 
out to be excluded if one assumes that the right-handed neutrino is of a third generation 
nature. In contrast, scenarios featuring a right-handed neutrino of one of the hrst two 
generations turn out to be excluded. In this way, we have demonstrated how accounting 
both for theoretical and experimental constraints has allowed us to almost exclude all 
possible setups in the left-right supersymmetric models under consideration. 

6 Conclusion 

We have studied a specihc class of supersymmetric models exhibiting a left-right gauge 
symmetry and investigated to which extent experimental and theoretical constraints restrict 
the viable regions of the parameter space. A particular property of this class of models 
that feature Higgs helds lying in the triplet representation of the SU{2)r group is that one 
of the corresponding doubly-charged Higgs bosons gets tachyonic when the scalar potential 
is minimized at the tree level. Even though it was known that loop corrections taking into 
account the Yukawa couplings of the right-handed neutrinos to the SU{2)e, Higgs bosons 
were modifying these conclusions, a complete one-loop calculation was still missing. We 
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have filled this gap and shown that in a large part of the parameter space, the complete 
one-loop calculation is necessary for reliable predictions. 

We have then studied to which extent the phenomenologically viable regions of the 
parameter space (where the doubly-charged Higgs boson is non-tachyonic), are constrained 
by experimental data and theoretical considerations concerning the global minimum of the 
one-loop-corrected scalar potential. We have found that the latter favours lower values of 
vji whereas the LHC gives tight lower bounds on this vev, in particular by exploring the 
Wr decays into jets. Moreover, the lightest doubly-charged Higgs boson is constrained to 
be well below 1 TeV for Yukawa couplings imposed to lie in the perturbative regime. The 
left-right supersymmetric setup that we have investigated is therefore close to exclusion, in 
particular if the future LHC searches for a Wr and a doubly-charged Higgs boson do not 
exhibit any signal within the next few years. 
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